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1 , BACKGROUND 

1.1 DEVELOPMENT OF THE DETECTION AND MAPPING (DAM) PACKAGE 

The beginnings of the DAM concept in the Spring of 1973 were pre- 
ceded by two events which occurred in the summer of 1972. The 
first was the launch of the Earth Resources Technology Satellite 
(ERTS) on July 25 r 1972. The satellite, then called ERTS-1 and 
subsequently redesignated Landsat-1, was the first of a series 
of Earth” looking satellites which provide relatively high resolu- 
tion multispectral scanner system (MSS) data to the general remote 
sensing community. 

The second event was the passage of Public Law 92-367 by the 
Congress of the United States and its approval by the President 
on August 8, 1972, two weeks after the launch of ERTS-1. This 
law charged the Secretary of the Army, acting through the Chief 
of Engineers,- "to . carry out a national program for the inspec- 
tion of impoundment dams. for the purpose of protecting human life 
and property." Implicit in this charge was the requirement to 
maintain an inventory of dams- at the national level , and to 
institute coordination or regulation, as required, of existing 
dam safety programs in the governments of all the states and 
territories. The initial effort at the inventory (1973-74) 
reported over 49 000 dams in the United States, each of which 

I 

impounded enough water to constitute a significant potential threat 
to the safety of areas downstream from the dam. The inventory 
data was compiled from many sources . The need for an economical 
and timely means for updating the inventory was recognized in 
the beginning. 

In December 197-2, the Texas Water Rights Commission (now the 
Texas Department of Water Resources) under contract with the 
Southwest Division of the U.S. Army Corps of Engineers, submitted 
through the Office of the Governor of Texas, a request for 
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assistance by the National Aeronautics and Space Administration 
(•NASA) Lyndon B. Johnson Space Center (JSC) , Earth Observations 
Division (EOD) in the development of a procedure for using 
Landsat data in detecting and locating water impoundments* In 
January 1973, representatives from the Office of the Chief of 
Engineers and NASA/JSC/EOD defined the following performance 
specifications for the requested inventory procedure. 

a. Landsat MSS to be the sole data source (no ground truth) 

b. Procedure to detect water bodies equal to or greater than 
ten surface areas to these criteria: 

o Correct identification equal to or greater than 90 percent 
o False alarms equal to or less than 10 percent 

c. Positional error less than 1000 feet 

d. Cost effective with conventional inventory techniques. 

The original research at NASA/JSC/EOD used existing general pur- 
pose multi spectral pattern recognition computer programs. These 
systems, designed for research and development applications, were 
not suitable for operational use. The situation led to the 
development of a special purpose hybrid procedure (utilizing 
both computer processing and human analysis) named the Detection 
and Mapping package. 

Periodic refinements have been made to the DAM package, which 
now consists of over ten main programs and more than 300 sub- 
routines together with associated graphic devices and procedural 
manuals. The unique software is directly useable on any Univac 
1100 series computer system under the standard EXEC- 8 operating 
system. The system characteristics of specific interest to the 
user are detailed in sections 2.4 and 3.0 of this manual. 
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In 1974, an early operational version of the DAM package was 
transferred to the State of Texas. The timing was such that this 
version was not used in the initial national inventory reported 
to the Congress by the Corps of Engineers . The DAM package was 
evaluated by two Corps of Engineers districts {Seattle and Nashville) 
who recommended its use for inventory update. These same two 
districts will now manage regional processing centers, providing 
computer-generated maps of surface water from the DAM package 
software for Corps of Engineers districts to use in updating the 
inventory of dams. 

1.2 NATIONAL PROGRAM OF INSPECTION OF DAMS (NPID) 

This program is directed by Section 2 of Public Law 92-367 which 
states in part; "All dams in the United States shall be inspected 
by the Secretary except (1) dams under the jurisdiction of the 
Bureau of Reclamation, the Tennessee Valley Authority, or the 
International Boundary and Water Commission; (2) dams which have 
been constructed pursuant to licenses issued under the authority 
of the Federal Power Act; (3) dams which have been inspected 
within the 12-month period immediately prior to the enactment 
of this act by a state agency and which the governor of such state 
requests be excluded from inspection; and (4) dams which the 
Secretary of the Army determines do not pose any threat to human 
life or property." Of all these conditions, the last one opens 
the subject of priorities. Priorities are needed not only because 
of the enormity of the task, but also because of the wide ranges 
of risks and hazards and the very real constraints of budgets. 

In a 1975 report to the Secretary of the Army, the Chief of 
Engineers advised that the cost of conducting an initial inspec- 
tion of a dam would vary from $5,000 to $10,000. The key element 
of all the coordination required with state agencies (in 1974, 

37 states had their own laws for dam safety programs) , and other 
federal agencies with jurisdiction over dams, is the dam inventory. 
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1.3 NATIONAL DAM INVENTORY 

A sample page of the inventory published in 1975 is shown in 
figure 1-1 of this manual. It is important for the analyst who 
will interpret .computer output from the DAM software to .be .aware 
of the essential elements of information required in the basic 
inventory. These include geographic, geometric, hazard, and 
ownership data. Because the analyst's functions are similar 
to those of aerial photointerpreta'tion , the analyst will be alert 
to other factors of interest in the area of dams. The most impor- 
tant of these are the nature of developments within the down- 
stream floodplain, and distances thereto. None of this informa- 
tion is contained in the computer output per se, but the analyst 
will inevitably ask questions about the area as soon as topo- 
graphic maps of an area are used. In any particular Corpse of 
Engineers district there are all kinds of area information in 
a variety of sources . It is reasonable to presume that the 
analyst, who is directly concerned only with the presence or 
absence of impounded water, will be involved with all of the data 
required to maintain the dam inventory. 

1.4 UTILITY OF THE DAM PACKAGE 

The maintenance of the National Dam Inventory (ideally in near 
real time) is clearly a formidable task. The special advantage 
of the DAM package is that it permits very rapid screening of 
large areas. The DAM computer processing discards all of the 
satellite data except detected surface water. The retained 
information is output by a computer line printer in map form, 
at any desired scale. It will be seen that these features will 
greatly accelerate the process of monitoring water impoundments 
on a nationwide basis . 
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2. AN OVERVIEW OF DAM COMPUTER PROCESSING 

The basic', purpose of the DAM package is the classification and 
mapping of surface water from remotely sensed MSS data acquired 
by NASA’s Landsat satellites. The information that ultimately 
emerges from the’ DAM package in the form of precisely registered, 
formatted, and interpreted maps begins as energy reflected from 
the surface of the ‘Earth. Small portions- of this energy are- 
detected and recorded high above the Earth by sensors aboard 
the orbiting Landsat satellites, as shown in figure 2-1. This 
data, now in digital or numerical . form, is transmitted -by the 

satellites to ground.-based receiving stations where, it is 

♦ ' 1 

recorded, retransmitted or carried to a centralized processing 
facility, computer-processed to apply numerous calibration 
and geometric corrections, reformatted, and distributed to users. 
The DAM package takes this data, in the form of magnetic tapes 
and photographic 'film transparencies, and further processes it 
to detect and precisely locate surface water and' to produce 
overlay maps depicting the surface water. These computer- 
generated maps .are, in turn, interpreted by hiiman analysts. 

This section presents a brief description of the chain of 
events that ultimately- leads to the creation of DAM package 
classification maps. 

2.1 LANDSAT MSS DATA 

Landsat is a system designed by NASA to provide for the repetitive 
acquisition of high' resolution MSS data of the Earth’s surface on 
a global basis. This overall system is illustrated in figure 
2-2. .To da'te, three Landsat satellites have been placed in 
Earth-orbit. Landsat-1 (now inoperative) was launched in July 
1972, Landsat-2 in' February 1975, and Landsat-3 in March 1978. 

Two sensor systems have been employed -to date on these Landsat 
satellites; a 4-channel (5-channel on Landsat-3) MSS and three- 
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Figure 2-2 Overall Landsat system . 
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(3 ON LANDSAT 1 AND 2 . 
2 ON LANDSAT-C) 


Figure 2-3.— Typical Landsat observatory. 
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camera (two-camera on Landsat-3) return beam vidicon (RBV) system. 
The elements of a typical Landsat observatory system are shown in 
figure 2-3. Since the DAM package utilizes only data from the MSS, 
the RBV system will be excluded from further consideration. 


The Landsat-1 

and -2 MSS* 

is a line scanning 

device that uses 

an 

oscillating mirror to continually scan perpendicular to the space 

craft direction of travel, 

as shown in figure 

2-4 . For each 


mirrbw sweep. 

six lines are scanned simultaneously in each of 

the 

following four 

spectral bands. 



MSS 

Landsat 


Wavelength 


channel 

band** 

Color 

(micrometers ) 


1 

4 

green 

0.5 - 0.6 


2 

5 

red 

0.6 - 0.7 


3 

6 

near 





infrared (IR) 

0.7 - 0.8 


4 

7 

near IR 

0.8 - 1.1 



The Landsat-1 and -2 satellites circle the Earth every 103 minutes 
in near-polar Sun-synchronous orbits at a nominal altitude of 
917 km (495 n. mi.) . They pass the Equator (southbound) at approxi 
mately 9:30 a.m. local time, complete approximately 14 orbits per 
day, and repeat the coverage of a given area once every 18 days. 
Figure 2-5 depicts some typical Landsat daylight ground traces. 

The Landsat MSS acquires data along a swath approximately 185 km 
(100 n. mi.) wide as the satellite circles the Earth. During the 
active time of one sweep of the MSS mirror, each of the four 
channels for each of six consecutive scan lines are sampled 


* Since the DAM package cannot presently accommodate Landsat-3 
data, only Landsat-1 and -2 will be considered in the 
remainder of this manual. 

**Bands 1, 2, and 3 are assigned to the RBV. 
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Figure 2-4.— Landsat MSS scanning arrangement. 
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Figure 2-5.— Typical Landsat daily daylight 
pass ground traces. 
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approximately 3240 times and the analog detector signals for each 
of these small square picture elements (pixels) are converted to 
digital form. Spacecraft forward motion provides for along-track 
progression of the scan lines. 

2.2 PREPROCESSING OF LANDSAT DATA 

The continuous -swath MSS data are transmitted to the ground where 
they are segmented into scenes comprised of exactly 2340 scan lines 
by approximately 3240 samples, which correspond to areas of 
approximately 100 by 100 nautical miles on the ground. Each 
scene is then further processed into various data products which 
are available to the public. The two such products utilized by 
the DAM package are listed below. 

a. System-corrected images — Each image covers an entire scene 
except for 42 lines at the top and 42 lines at the bottom. 

The scenes are radiometrically calibrated and nominally 
corrected for scanner geometry, spacecraft attitude/altitude 
and Earth's rotation based on satellite tracking and tele- 
metry information. One scene is comprised of approximately 
7 500 000 individual pixels, each representing, (by a 
density or shade of grey) the radiance reflected, for the 
particular spectral band from a 79- by 57-meter area of the 
Earth's surface. Figure 2-6 shows a typical system- corrected 
Landsat MSS image of MSS channel 4 (Landsat band 7) . 

b. System-corrected computer-compatible tapes (CCT's) — Digital 
MSS data for each 100- by 100-nautical mile Landsat scene 
are divided into four 25-raile wide strips. The data are 
radiometrically calibrated but are not corrected for scanner 
geometry, platform attitude variations, orbital altitude, or 
other perturbations. The data for the four str-ips of a scene 
are recorded on either one, two, or four reels of magnetic 
computer tape, depending on the format and recording density - 
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2-6.— Typical system-corrected 
image of Landsat scene. 
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2.3 DAM PACKAGE DATA PROCESSING 

The DAM software utilizes the surface radiance of Earth features, 
as recorded by the four Landsat MSS channels to classify or 
detect surface water. The digital data used to perform this 
classification are geometrically uncorrected and in a dynamic 
scanner-oriented coordinate system. After classification these 
data are corrected and transformed into an Earth-based coordinate 
system. The correction and transformation process, termed 
registration, utilizes parameters derived from a network of 
"control points", identified and measured on both maps and the 
MSS scene, to bring the classified data into registration 
with prescribed base maps. 

The generalized flow of the data processing within the DAM package 
including control network establishment, surface water detection, 
and generation of surface water maps, is summarized in figure 2-7. 

2.4 DAM PACKAGE OUTPUT DATA FORMAT 

The standard DAM package output maps are produced by a conven- 
tional, 132-column, computer-driven line printer. Although the 
DAM software can output correctly scaled maps with a paper feed 
of six lines per inch (standard commercial) , better map 
resolution is achieved when the printer feed is switched to 
eight lines per inch. This section discusses the general aspects 
of the output line-printer thematic maps. Specific elements 
of the maps are detailed in section 3. 

2.4.1 LINE PRINTER MAPPING GEOMETRY 

The line printer outputs one horizontal line of discrete characters 
at a time, and then advances the paper vertically before the next 
line is printed. At eight lines per inch and ten columns per 
inch, each print character effectively covers an area 0.125 inches 
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2-7.— Generalized DAM package data flow. 
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Figure 2-7.— Continued. 
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Figure 2-7.— Concluded. 
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high by 0.10 inches wide. Now recall that the Landsat MSS 
acquires data by sampling the reflected radiation for discrete 
pixels within scan lines, and that orbital motion advances the 
satellite to the south before subsequent scan lines are sampled. 

It is a relatively simple matter to assign each pixel to one 
print character. In this case, there is a direct correspondence 
between print lines and scan lines, and also between printer 
columns and scanner samples. 

Several factors complicate this situation. First, the proportions 
of height to width for print characters and scanner pixels are not 
the same. Second, scan line and scanner sample are not at right 
angles to each other due to the skewing effects of Earth rotation 
and yaw. Third, the area on the ground (and hence, the number of 
pixels) represented by a constant-size print character varies 
with the scale of the map. 

Figure 2-8 illustrates the way in which the registration routines 
within the DAM software handle these complications when generating 
line-printer maps at a scale of 1:24 000. The rectangular grid 
outlining the print characters is sketched in solid lines and the 
skewed grid outlining the scanner pixels is sketched in dashed 
lines. The center of each pixel is marked with a dot. Note 
the number of pixel centers that lie within each print character: 
some characters have no pixel centers, most have one, and some 
have two. Clearly, this results from the proportion differences, 
skewing, and scaling described in the previous paragraph. By 
holding figure 2-8 up to the light you will see the symbols that 
the registration software would assign to these print characters 
if all the pixels were classified as water. 
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Figure 2-8.— Locations of Landsat pixel centers 
relative to line printer character locations 
for 1:24 000-scale computer-generated maps. 
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The meaning of these symbols is as follows; 


Symbol 


Meaning 

0 pixel centers (no data) 


1 


1 pixel center 


2 


2 pixel centers 


Because (at 1:24 000) a print character is shorter than a pixel, 
periodically a line of print characters will occur for which there 
are no corresponding pixel centers. Any of these print characters 
below pixels classified as water will then be printed as a colon 
(;) and any not below water pixels will be printed as a blank. 
Conversely, (at 1:24 000) a print character is wider than a 
pixel, resulting in the periodic occurrence of print characters 
with two pixel centers, symbolized by an arabic numeral 2. Note 
that these 2 ' s occur on the diagonal as a result of the previously 
described skew in the placement of pixels. 

The net result of these registration techniques is that pixels 
classified as water are never duplicated (a no-data symbol is 
printed instead) nor are they lost (an arabic numeral 2 is 
printed instead). Hence, an exact count of pixels classified as 
water can always be determined from the registered output maps . 

At scales other than 1:24 000, the same techniques are used. 

Each print character contains a symbol indicating the number of 
pixels classified as water whose centers lie within the character. 

I 

An actual line printer example of these effects is shown in 
figure 2-9. 
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Figure 2-9.— Large water body at 1:24 000 showing 
line of NO DATA colons and periodic "2s". 
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2.4.2 MAP ORIENTATION AND TICK MARKS 

Efficient computer implementation of the registration techniques 

f 

explained in the previous section dictates that print lines on 
the line printer output maps must be parallel to scan lines in 
the MSS data. For this reason / individual line printer maps are 
oriented to the local orbital direction of the Landsat satellite. 
Because the Landsat orbital plane is inclined at 9.3 degrees to 
the Earth's axis (to achieve Sun-synchronous coverage) the MSS 
scan lines deviate by 9 . 3 degrees from true east-west at the 
Equator. This angle increases at higher latitudes. Figure 2-10 
shows the typical orientation of DAM package computer-generated 
map relative to a standard USGS quadrangle map. 

Not shown in the figure are the marginal data and the registration 
tick marks printed on the classification map. The DAM software 
computes the line printer location of all four corners of the base 
map. These are marked with asterisks (*) and designated primary 
tick marks. Secondary tick marks, marked by a plus sign (+) , 
are plotted at even grid intervals on the map inside the four 
corners. The standard tick interval for the 7.5-minute USGS quad- 
rangle map sheet is 2.5 minutes in both latitude and longitude. 
These tick marks assist the analyst in rotating and shifting a 
computer-generated map to match the corresponding conventional map. 

2.4.3 MAP COVERAGE 

The DAM mapping routines automatically calculate the number of 
line printer units that will be required to accommodate a particu- 
lar map at the specified size and scale. For the standard USGS 
quadrangle map sheets, the size in the vertical (north-south) 
direction remains constant at 7.5 and 15.0 nautical miles (1852 
meters/nautical mile) for the 7.5-minute and 15-minute quadrangle 
sheets, respectively, at their respective scales. Map length 
in the north-south direction is generally in the direction of the 
printer paper feed, except that as the latitude increases, the 
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Figure 2-10.— Orientation of DAM package map data 
with standard USGS map sheet. 
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Landsat data orientation rotates with the changing spacecraft 
heading, in a westerly direction. For example, at the western 
border between the United States and Canada at 49° N latitude, 
the angle between the Landsat data and meridan lines on the map 
increases to nearly 14 degrees. On the other hand, the base map 
width has decreased in proportion to the cosine of the latitude. 
The two effects counteract each other. It is enough for the 
analyst to know how to determine how many computer printouts 
should have been prepared for each map overlay. This is explained 
in section 3.1. 

2.4.4 LANDSAT DATA OVERLAP 

The Landsat orbit progression and the field of view of its sensors 
have been designed so that, in the worst case at the Equator, 
the nominal data sidelap between longitudinally adjacent frames is 
7.5 nautical miles. At all other latitudes sidelap is greater. 

It increases in inverse proportion to the cosine of the latitude. 
One minute of longitude is equivalent to one nautical mile at 
the Equator, therefore, the need to merge data from two or more 
scenes is reduced since one scene or another will usually contain 
all of a 7.5-minute quadrangle map sheet. Larger maps could 
present a problem, but it is unlikely the 15-minute quadrangle, 
for example, can do so in temperate latitudes and certainly not 
in extreme latitudes. As the satellite data overlap increases, 
the longitudinal distance of the ground covered by the map 
decreases . 
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3. PREPARATION OF COMPUTER-GENERATED MAPS 


A recurring step throughout the actual interpretation and 
analysis of DAM package water classification maps will be the 
comparison of that information to the related topographic map. 
The simplest way to make these comparisons is by overlaying and 
registering the two maps, and plotting any other information on 
the overlay (or the map) . When one line printer printout makes 
a complete map, this one map unit is easily registered to the 
base map for which it was constructed. In the more usual case, 
at least two separate computer printouts must be assembled to 
show water classification results over the entire base map area. 
When this occurs, the second (right most) unit must be carefully 
trimmed so that it can be joined in exact alignment with the 
first unit, with no loss of essential marginal data. This is 
called a mosaic. The mosaic may then be registered to the 
base map. The construction of a three-unit mosaic is not 
described in this manual. It will become obvious that the 
three-unit map is joined by simply repeating the process used 
to construct the two-unit map. 


3.1 UNPACKING THE DATA 

When data covering one Landsat scene arrives from the computer 
center it will be in the form of several stacks of computer print- 
outs. The cover or header sheet of each DISPLAY /MAP PRINT FILE 
will be similar to the one shown at the top of figure 3-1. The 
top line of block symbols are computing system identifier codes. 

The next line identifies DAM PRTCLASS followed by a version number, 
(DAM PRTCLASS is the computer program which produces the registered 
maps on the line printer) . Then follows the date and time of the 
processing; and finally, the file identifying numbers. The left 
block n\imber identifies the number of the first file (0) , while 






3-1 


VERSION 7802 


JSC-13970 




LCLCl 

sss 

L 

tsa 

X 

X 



ClCLL 

SSS 

L 

888 

X 

X 





L 

s s 

L 

t e 

X 

X 


/ 

L 

a s 

L 

e b 

X 

X 





£ 

s 

L 

8 a 

X X 


/ 

fa 

s 

L 

b 8 

X 

X 





tCL 

sss 

L 

tBS 

X 


/ 

sss 

L 

8 86 

X 






i 

s 

L 

6 8 

X X 


/ 

c 

s 

L 

b b 

X 

X 





£ 

s s 

L 

6 8 

X 

X 


/ 

c 

s s 

L 

a a 

X 

X 





LCLCL 

sss 

LLLLL 

fesa 

X 

X 



ClCCC 

sss 

LLLLL 

8b8 

X 

X 



LLDG 


AAA 



Phpp 

RRRR 

T IT II 

ccc 

L 

AAA 

SSS 

SSS 


1 

77777 

888 

L 

c 

A A 

MH HM 


p p 

R 

R 

T 

C 1 

1 

A A 

S S 

s S 



7 

8 8 

L 

c 

A A 

»* •» N 


p p 

R 

R 

I 

C 

L 

A A 

S 

s 



7 

8 8 

U 

0 

AAA AA 

»• •• M 


pppp 

RRRR 

I 

C 

L 


>ss 

sss 



7 

888 

L 

0 

A A 

N m 


p 

R R 

T 

C 

L 

A A 

s 

s 



7 

6 8 

L 

0 

A A 

p p 


p 

R 

R 

T 

C L 

L 

A A 

s s 

s > 



7 

8 8 

LGOl 


A A 

p p 


p 

R 

R 

T 

CCC 

LlLLL 

A A 

:»SS 

sss 


< 

7 

888 


00 


222 I 
2 2 
2 

222 

2 

22222 I 


S.0 

0 0 
C 0 
0 9 

0 0 


•I 

<1 


1 

11 


1 

111 


CO 


77777 

7 

7 

7 

7 

7 

7 


6 a 

c6tt 
fe 8 

» a 

«8b 


-0 

L 0 
0 3 
5 9 
o 0 


333 
3 3 
3 

33 


222 
2 2 
2 

222 

2 

22222 


333 
3 3 
3 

33 

3 

3 3 
333 


«.lSf LA V /MAP OUlPOl 


UO Cu 

CGOOC^OC 


LOC 

CbO 

CO 

00 

CO 

CO 

cuo 

COG 


0 uC 

coo 

§s 

CO 

CJ 

LOO 

CLC 


CGC0C9CG 
COCO 


••• iHlb ••• 

11 

All 

1111 

11 

li 

11 

11 

ii 

11 

11 

11 

11 

iimi 

iiiiii 

• •• 71tL •• • 


LASI 

333333 
3333333333 
333 333 

33 33 

33 

333 

333 

n 

ill 


33 
133 
3333333333 
333333 


r ILL 


«•« 


C!»LrV /LSL6X 
ICb^C-lto3G3 / X0C976 


uA>- PPTCLASS «7b*'«. I 
/ USAIO SCAlTLc. y lHS 


Z*t/*C/7L u3:23 HlE •LAH^RI-I 


• NO OrCK • 


?7?‘»??7 ????????????????????? ???????? 


LC YOU H17C : 

1 COPT OF BASK PRIM flLf FOR IwIS MUK??? 

1 COPUS or ALL A GISPLAV/HAP PMIM FILlS fop This PPCSkAh CxtCuTlOL??? 
??7??7???r? 7?? ????????????????????? Y 




Figure 3-1.— Typical DAM package 
DISPLAY/MAP PRINT FILE cover sheet. 
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the right hand block identifies the sequence number of the last 
file (in this example, 3). The large block number centered at 
the bottom of the page (middle of fig. 3-1) is the number of this 
particular file (1) . Thus, the file in this example is the second 
of a total of four separate files (0 through 3) which make up this 
data set. The total number can be more or less than four. This 
will depend on how line printer operations are conducted at the 
computing center serving you. 

3.2 COLLATING THE COMPUTER MAPS 

Place the four numbered files side by side, in ascending order 
from left to right, in the manner indicated in the top line of 
figure 3-2. Remove the header sheet (and its backup sheet which 
contains a summary of the number of files and number of copies 
each) from each file, which will expose the top page of a map 
unit in each file. 

The UNIT ZERO of the lowest numbered map (normally 001) , will be 
on top of the first file at the left. Figure 3-3 is an example 
of this page. UNIT ZERO (if present) normally consists of three 
or four fanfold pages as depicted schematically in figure 3-2. 

These pages contain formatted annotation data, directly analogous 
to the marginal information on a conventional map. The UNIT 
ZERO should be torn off from the rest of that file at the end of 
its last page (fold the pages back into their original orienta- 
tion in the file) . The graphic part of the map associated with 
this UNIT ZERO consists of one or more subsequent units (UNIT 
ONE, UNIT TWO, etc.). 

The UNIT ONE for this map will be on the top of the file immediately 
to the right of the file which contained UNIT ZERO. The UNIT ONE 
may now be detached from its file (split the perforations between 
the last sheet of this unit and the header sheet of the next unit) . 
Refold this unit and place it under the UNIT ZERO for this map 
which was previously detached. 
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Figure 3-2.— Typical sequence of maps and units 
within DISPLAY/MAP PRINT FILES. 
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Repeat the above process to find, separate, refold, and collate 
the UNIT TWO (and THREE, if it exists - for standard USGS quad- 
rangle map sheets, there is no possibility of a UNIT FOUR) . 

Once all units of a given computer map have been collated, continue 
the process by finding the UNIT ZERO tif present) of the next 
higher numbered computer map (normally 002) . According to the 
pattern, this new UNIT ZERO should appear on the top of the file 
immediately to the right of the file which contained the last unit 
of map 001 » Figure 3-2 illustrates how the nximber of units varies 
from map to map. Each computer map will always have a UNIT ONE 
and may or may not have a UNIT ZERO. The printer page length of 
the various units varies from map to map. This separating, 
refolding, and refiling process may now continue until all of the 
original four files have been redistributed into however many 
integral computer maps they contained. 

3.2.1 COMPUTER MAP BOUNDARIES 

The completeness of a DAM computer map can be verified by examining 
the boundary frcimes of UNIT ONE (and TWO and THREE if present) . 

Each unit map area is bounded on all four sides by a single line 
(column) of colons. (The colon is used throughout as the symbol 
for no data.) These colons correspond to the neat lines on a 
conventional map. Immediately above and below the colons at the 
top and bottom of the computer map there are four-digit printer 
column numbers. These increase from left to right. Immediately 
outside and adjacent to the colons bounding the sides of the map 
area there is another series of four-digit numbers. These are 
printer line numbers. These increase from top to bottom. 

If a map UNIT ONE has line numbers printed down both sides, then 
there will be no corresponding UNIT TWO. If UNIT ONE has colons 
only, without numbers down the right side boundary, then there 
is a continuation of the map data on a UNIT TWO. If there 
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should also be an associated UNIT THREE, the UNIT TWO will have 
only colons along the right margin, and so on. These character- 
istics provide the analyst with an independent check on the 
completeness of any particular DAM map unit set. Figure 3-2 
illustrates schematically the side boundaries of a two-unit map. 

3.3 CONTENTS OF THE UNIT ZERO 

This unit contains the computer map legend, including information 
on source, reliability, and scale. Because UNIT ZERO information 
is not needed for all computer maps, UNIT ZERO sheets may be 
selectively suppressed by the regional processing center. 

3.3.1 PAGE ONE 

Refer to figure 3-3. In the top line, the numbers following DAM 
PRTCLASS give the software version, processing date, and time. 

The last triplet of numbers identify: the computer map num- 

ber within this program execution - the unit number within this 
computer map - and the page number within this map unit. The 
second line gives accounting, routing, and quality control in- 
formation. The third line normally gives the map name, series, 
and state. 

The large block ID numbers repeat the time, map number, and 
unit number, in a form easily read from a distance. 

Scene date is the date on which the scanner data was recorded 
by the satellite. Map origin for a standard map (conventional 
or computer-generated) is normally the latitude and longitude 
of its southeast corner. 

3.3.2 PAGE TWO 

This page presents a legend of the symbols used on this map (see 
fig. 3-4). The NO DATA (:), PRIMARY TICK (*) , and SECONDARY 
TICK (+) symbols are always used for these purposes, and for no 
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Figure 3-3.— Typical UNIT ZERO page 1. 
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Figure 3-4.— Typical UNIT ZERO page 
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other. The remaining symbols describe the ntxmber of water- 
classified Landsat MSS pixels whose centers lie within the line 
printer character displaying the symbol. 


3.3.3 PAGE THREE 

Figure 3-5 illustrates this page. It contains part of a table 
giving coordinate locations for the tick symbols printed on 
subsequent units of the map. The exeunple shown is for a 15- 
minute, 1:62 500-scale map. 


Tick tables are continued on additional paper, if required. 




3.4 CONTENTS OF UNIT ONE (AND HIGHER) 

The marginal information presented at the top of UNIT ONE is 
identical to the presentation on page one of the UNIT ZERO with 
the single exception of the unit number at the end of the very 
top line, and in the last of the large block ID numbers. The 
rest of the printout is all map data, surrounded by the four 
line printer neat lines as described in section 3.2,1. Figures 
3-6 and 3-7 show the units associated with map 026. There are 
three special rules concerning the appearance of tick symbols. 

a. Primary tick symbols are always printed on the computer 
map, even at the expense of suppressing (eliminating) water 
data which would have been present in the character 
occupied by the asterisk . 

b. Secondary tick symbols are printed at all proper locations 
on the computer map (these can be confirmed by noting the 
line printer coordinates in the tick table, UNIT ZERO) 
provided they do not displace either a water symbol or a 
primary tick symbol. Water symbols take precedence over 
secondary tick symbols . Any secondary tick symbol which has 
been suppressed on the computer map due to conflict with a 
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Figure 3-5.— Typical UNIT ZERO page 3. 
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Figure 3-6.— Top of typical UNIT ONE. 
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Figure 3-7.— Top of typical UNIT TWO, 
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water symbol may be manually plotted, if necessary, from its 
coordinates as given in the tick table. 

c. Both primary and secondary ticks appear in place of no data 
symbols. Further, the no data symbols immediately adjacent 
to the tick symbol on that printer line are also suppressed, 
providing a "halo" which assists in locating the tick mark. 
These special situations appear in the NO DATA area at the 
top of the computer map shown in figure 3-6. 

A final comment on the computer-generated maps — because they are 
constructed at an oblique angle with the neat lines of the re- 
lated base map, they will always subtend a larger dimension, 
both vertically and horizontally, than the related base map. 

3.5 INDEXING AND SCREENING 

Once the collating process is completed, the data is organized 
into sets of computer-map sheets ("units"), ready for assembly 
into the individual complete computer maps (see section 3.6). 

At higher latitudes and smaller scales, the complete map will 
probably consist of just one unit. At lower latitudes and 
larger scales, the complete computer map will be composed of 
several units. A full Landsat scene will generate between 200 
and 300 maps at various scales. Fewer maps will be produced 
for scences which, in part, are cloud-covered or lie outside 
the area to be mapped. The largest number of computer maps 
will be those which cover the scene area at the largest scale, 
that of the USGS quadrangle map (1:24 000). 

The lowest numbered computer map(s) for a Landsat scene normally 
will cover the entire usable portion of the scene at a scale 
of 1:250 000, with tick symbols every 7.5 minutes of latitude 
and longitude. These 1:250 000-scale maps provide an over- 
view of the classification results for the scene. They are 
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useful as an indexing tool for the larger scale computer maps of 
the same area, and for identifying areas of potential diffi- 
culties because of haze or clouds, by reference to an original 
Landsat data scene. For partial Landsat scenes, "overview" 
maps may be omitted. 

The 1:250 000-scale computer-generated map(s) should be assembled 
first, before any other maps for this scene, using the techniques 
given in section 3.6. Once this has been done, use colored 
pencils and straightedge to outline the boundaries of 15-minute 
(blue) and 7.5-minute (red) quadrangle maps by drawing the 
meridians and parallels connecting all the tick symbols. The 
latitudes (longitudes) of these lines should be determined from 
the tick table in UNIT ZERO, and then annotated on the map. Once 
this is completed, the names of the published quadrangle maps 
can be determined from the standard USGS 1:1 000 000 index map, 
and then annotated on the 1:250 000 computer-generated map. 

This annotated map will now be used along with other considera- 
tions, to determine which individual large-scale computer maps 
need to be assembled and the order in which they should be 
interpreted. 

The first consideration is the areas (s) assigned for interpreta- 
tion and analysis and perhaps priorities. This is an operational 
matter which will vary between districts. The technical considera- 
tions are the availability of maps and the location of areas 
which will require special attention because of spotty cloud 
cover, or which are unusable because of extensive haze or cloud 
cover. Normally the latter will have been eliminated by the 
regional processing center. The regional processing center will 
process only the most cloud-free scene available. An overlay show- 
ing the published map coverage in the area will be furnished by 
the Regional Processing Center, Cloudy or hazy areas will usually 
be indicated on this overlay. The analyst should examine the over- 
lay as a check, and then annotate the 1:250 000-scale computer map. 
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The areas of affected maps to be interpreted with caution are now 
well pre-identif ied. In cases where proper evaluation cannot be 

i 

made, the regional processing center will provide alternate 
coverage on request. 

The larger scale computer maps to be assembled can now be selected. 
Computer maps are selected at scales which match the large scale 
topographic maps available. This applies to all computer maps, 
including those which indicate no detected water at all. After 
assembly, each computer map will be overlaid on its related 
topographic map to check for mis-located or apparently missing 
dams, as well as new impoundments. Use of the largest map 
scales available assures the best accuracy possible of interpre- 
tation and coordinate measurements. 

The indexing and screening has now been completed. 

3.6 TRIMMING AND FITTING THE MOSAIC 

The joining of two contiguous aerial photographs into one is 
called mosaicking. The same process is followed in assembling 
UNIT ONE and UNIT TWO into a complete computer map. It is made 
easier by the fact that the line printer character geometry 
automatically provides a grid which is used as a guide for 
trimming the excess paper away from the left edge of UNIT TWO 
and for positioning it exactly tangent to the right edge of 
UNIT ONE. Figure 3—8 shows where to cut off the left margin of 
UNIT TWO. The cut is made vertically from the very top of the 
unit to the bottom, immediately to the right of the column of 
colon symbols (:) which signify the left edge of the mapped data. 
The column to the right of the colons contains map data. When 
it is blank, the column can still be located by noting the 
print column scales at the top and bottom of the map area. 
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Figure 3-8.— Illustration of location 
of trim line for UNIT TWO. 
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The UNIT TWO sheet may then be fitted over the right margin of 
the UNIT ONE sheet. Long heavy straightedges or map weights 
should be placed along the entire length of the join line, on both 
sides, to ensure that the map units lie perfectly flat. Once both 
units are aligned so that the print column scales at the top and 
bottom line up exactly, both vertically and horizontally, they 
should be taped together securely. 

i 

It will be noted that trimming the UNIT TWO sheet rather than the 
UNIT ONE produces an overlap joint with no printing on the under- 
laying map edge (other than a single column of colons) . If the 
UNIT ONE were trimmed and fitted over the UNIT TWO, the UNIT TWO 
left margin line numbers would show through on a light table, 
making registration and interpretation more difficult in the 
region of the mosaic seam. 

The marginal data provided at the top, bottom, and sides of the 
UNIT ONE/UNIT TWO mosaic are important and must not be trimmed 
off. Similarly, the legend and tick table in UNIT ZERO are 
important and should be kept with the associated mosaic. 
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Before examining the details of interpretation and analysis of 
DAM computer maps, it is useful to recall the criteria for 
reporting in the national inventory. These are: all dams 

25 feet or higher, which impound at least 15 acre-feet of water, 
and all dams which impound 50 acre-feet or more, which are at 
least 6 feet high. Obviously, satellite data alone will never 
provide measurements of this precision. The whole purpose of 
the DAM package is to provide a method of rapidly screening 
very large areas for surface water in general , and impounded 
water in particular. Discrimination between impounded and non- 
impounded surface water, location of the impounding structures, 
and the reduction of misclassif ications caused by gross shadows 
are the key interpretation services provided by the analyst. 

4.1 PIXEL SIZE AND CONTIGUITY 

It will be recalled that at a scale of 1:24 000, one Landsat 
pixel is approximately equal in area to one line printer character. 
The area on the ground uniquely recorded by one Landsat pixel 
covers approximately 1.1 acres. The DAM package classification 
processing selects "pure" water pixels, which means that the 
pixels classified as water are entirely water, and contain no 
significant soil, rock, or vegetation on the surface. Pixels 
which contain mixtures of water and vegetation (typical of 
marshes and of lakes choked with hydrilla) , or water, shoreline 
and land combined, are classified as non-water. For these 
reasons, pixels classified as water will generally be surrounded 
on all sides by edge pixels not classified as water. 


\ 
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The number of edge pixels and the proportions of water to non- 
water in each edge pixel will vary with the size and shape of 
the water body as well as the degree of alignment between its 
shoreline and the grid of scanner pixels. As a consequence, 
the surface area of a water body cannot be determined by simply 



counting the number of contiguous pixels classified as water 
and multiplying by 1.1 (or by any other number) . For example, 
three contiguous pixels classified as water generally correspond 
to a surface area of approximately ten acres , but the actual 
area of any individual 3-pixel water body may vary from four 
acres to twenty acres. Likewise a single pixel classified as 
water corresponds to a surface area between two and ten acres. 


Nevertheless, the number of contiguous water pixels detected is 
a useful, if crude, indication of the areal extent of the actual 
water body. Contiguous pixels are defined as those which are 
adjacent at a side or corner. 



Pixel contiguity can always be determined from the 1:24 000 scale 
computer maps. The line printer symbols will be adjacent at a 
side or corner if and only if the corresponding pixels are 
adjacent. In this determination, two water symbols represent 
contiguous pixels even if separated by a mutually contiguous 
NO DATA (:) symbol. 

Pixel contiguity cannot always be determined from computer maps 
at 1:62 500 and smaller scales. Figure 4-1 illustrates some 
problems in attempting to determine pixel contiguity from dif- 
ferent scale computer maps. At "A”, two symbolizations of 
three contiguous pixels at 1:24 000-scale are shown, and 
several different ways these can appear in the related 
1:62 500 computer map. 
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At "B”, two symbols (and their corresponding pixels) are con- 
tiguous, but the third (and its corresponding pixel) is not. In 
the 1;62 500 computer map, two different presentations are shown 
of the same data (their difference lies in where the 1;62 500 map 
character falls with the actual Landsat pixel centers) , Without 
the information contained in the 1:24 000 map, the analyst would, 
looking at either of the 1:62 500 presentations, erroneously 
conclude that all three pixels were contiguous. 

Part "C” of figure 4-1 illustrates one other special situation. 
The 1:24 000 computer maps show water areas which include 
NO DATA colons inserted by the computer program to preserve 
the map scale (see section 2.4.1). Water symbols which are 
above and below a colon, or diagonally opposite across a 
colon, ^ represent contiguous water pixels. 

Of course, the final determination of whether neighboring non- 
contiguous pixels belong to two separate water bodies or to a 
single very narrow water body cannot be made by any computer 
but depends on human judgement in interpreting both the computer 
maps derived from satellite data and the elevation contour lines 
on topographic maps. 
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4.2 PRELIMINARY EXAMINATION OF THE DATA 

In the two following examples, USGS 15-minute quadrangle maps 
were available, but not 7.5-minute quadrangles. The DAM software 
printed out a single 1:62 500 map and the four subordinate 
1:24 000 maps. Prior to beginning interpretation of the data, 
the analyst should study the topographic map to familiarize him- 
self with the general nature of the terrain. He should also note 
the year in which the map was compiled and field-checked, and 
compare this with the date of the Landsat data. To commence in- 
terpretation, the analyst places the topographic map and the 
corresponding computer-generated map on a light table. 

4.3 REGISTRATION TO THE BASE MAP 

The computer generated map and the conventional map, overlayed on 
the light table, are shifted and rotated so that the computer- 
produced tick symbols match the tick marks on the conventional 
map. 

It is important to note that standard 7.5- and 15-minute USGS 
quadrangle maps occasionally are extended (usually in the vicinity 
of international boundaries) to cover areas outside of their nor- 
mal corner tick marks. In these special cases some corners of 
the published maps will not correspond to tick marks; the analyst 
must locate the actual ticks using the annotations in the 
quadrangle map margins. 

Which map should be on top? If there is much apparent surface 
water, the analyst will probably prefer to have the topographic 
map on top of the computer-generated map. The elevation contour 
lines describe the terrain relief, a major consideration in 
ruling on both terrain shadows and whether or not a water body 
is natural or impounded. Since the contour lines are printed in 
a moderate brown color in the standard four color map, these 
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must be viewed directly if detailed study is required. However, 
if the interpretation can be made with the computer map on top, 
this arrangement is more convenient for noting conclusions 
directly on the computer map. 

4.4 ANCILLARY MAPS 

If published 1:24 000 topographic maps are not available, the 
1:24 000 computer-generated maps need not be assembled. How- 
ever, their units must be available to the analyst since they 
may be needed to examine apparent water impoundments with better 
resolution. The 1:24 000-scale computer maps which cover the 
four 7 . 5— minute quadrants of a 15— minute map immediately follow 
the 15-minute computer map. Thus, if the 1:62 500 map is number 
015, then map numbers 016-019, inclusive, are the 1:24 000 maps 
for its four 7.5-minute quadrants. These four 7.5-minute 
computer maps will all bear the same name (that of the 15-minute 
quadrangle) , but can be distinguished by the latitude and 
longitude of their origins (southeast corners) . 

4.5 REJECTION OF SPURIOUS WATER PIXELS 

As discussed in section 3.5, inspection of 1:1 000 000 Landsat 
image or overlay furnished by the regional processing center will 
quickly show whether the areas to be analyzed are cloud free. If 
they are clouded, the regional processing center should be asked 
for alternate data. However, if the clouds are scattered and 
relatively small the data is marginal and may be useable with 
caution of two kinds - interpretation of the potential for a water 
impoundment in the terrain obscured by the cloud itself, and 
identification of apparently detected surface water which is 
actually the cloud's shadow. An apparent water body located on 
the sunny slopes of a mountain without a source watershed may well 
be a cloud shadow. The cloud itself will be located along the 
back-azimuth of the sun azimuth, and will be about the same size as 
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the shadow. The distance may be anywhere from about one to 
three miles (assuming typical cumulous cloud altitudes of 
10,000 feet, sun angles from 30° to 60°) . If terrain factors 
suggest possible impoundments in the vicinity of clouds, the 
data should be rejected. 

In the downtown areas of major cities, other spurious indications 
of water from the satellite data can arise: large groupings of 

tall buildings will cast shadows on the streets and parking areas 
to their northwest. If these are extensive enough and dark 
enough, some of the pixels falling within them will likely be 
indistinguishable to the Landsat MSS from clear, deep water. 

Most such misclassif ications (as locating a lake in the midst of 
the New York financial district) can be readily screened out by 
the analyst. However, some caution is necessary: lakes within 

urban parks will be properly classified as water (if large enough) , 
and are likely to be man-made (and thus, reportable). Likewise, 
cooling ponds, settling basins, sludge ponds, acid pits, flooded 
quarries, abandoned contour strip mines, large flat building roofs 
with standing water, etc. will be correctly identified as surface 
water by the computer processing (if large enough) , and some of 
these may be reportable. 

Transmission problems with data from the satellite may occasion- 
ally result in straight horizontal or vertical lines of apparent 
water running across an entire scene, particularly along its 
rightmost edge. These spurious "water bodies" should be ignored. 

We now move to an analysis of spurious water detection caused by 
terrain shadows. Figure 4-2 is the northeast corner of the Van 
Zandt 15-minute quadrangle map in Whatcom County, Washington/ 
on the western side of the Cascade Mountains. The terrain is 
extremely rugged and the annual rainfall is substantial. The 
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Figure 4-2.— Northeast corner of Van Zandt 
(Washington) 15-minute USGS quadrangle map. 
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Figure 4-3.— DAM package computer-generated map 
corresponding to portion of USGS 15-minute 
quadrangle -map shown in figure 4-2. 
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Figure 4-4.— DAM package computer-generated map of 
figure 4-3 superimposed on portion of Van Zandt 
15-minute USGS quadrangle map shown in figure 
4-2. 
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developed areas are confined to the narrow valleys of the Nooksack 
River and its tributaries. The Canadian border is the top neat 
line of this map.* Figure 4-3 is the corresponding part of the 
related computer-generated map. the marginal information on this 
map indicates that the data was obtained by Landsat-1 on September 
18, 1973, and that at the time the data was acquired, the Sun 
elevation was 39 degrees and the Sun azimuth was 149 degrees. 
Therefore, some northwest slopes in the steep terrain of Van 
Zandt area may be shadowed. Examination of figure 4-4 proves 
this to be so. The water pixels confirmed on the map are (1) 
a natural lake, Silver Lake; and (2) the meandering portions of 
the Nooksack River, big enough to contain "pure” water pixels. 

All of the other "water" pixels which have been circled are clear- 
ly not water. They are all located in the steep areas of north- 
west slopes, or in deep draws or ravines. In 1973, then, there 
were no impoundments in this area where people would be extremely 
vulnerable to an unsafe dam. This data set illustrates the 
analyst's critical role: to add intelligence to the DAM package 
which the computer cannot. 

4.6 SELECTION OF QUALIFYING WATER BODIES 

To further illustrate this process of interpreting water bodies, 
more data is presented from the same Landsat scene. The area 
selected is 15 nautical miles south of the Van Zandt quadrangle. 
Figure 4-5 shows the northwest area of the Clear Lake 15-minute 
quadrangle. The principal features of this region are, again, 
the Cascade Mountains, but also more lowlands. The developed 
community of Sedro-Woolley is partially included in the quadrangle 


In section 4.3 it was suggested that not all USGS quadrangle 
sheet corners would be at exact multiples of 7.5 minutes or 
15 minutes. This sheet is one of the odd ones. It extends to 
the north about 200 meters beyond the 49th Parallel to the 
actual United States-Canada border in this area. 
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Figure 4-5.— Northwest portion of Clear 
(Washington) 15-minute USGS quadrangle 
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Figure 4-6.— DAM package computer-generated map 
corresponding to portion of USGS 15-minute 
quadrangle map shown in figure 4-5. 
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map. Other major features are the Skagit River and several 
natural lakes. Because this data is from the same Landsat scene 
onal a few miles south on the same Landsat orbit, the Sun angles 
will be identical. (Landsat overflies an entire 100-nautical 
mile scene in about 26 seconds.) 

Figure 4-6 is the corresponding portion of the computer-generated 
map. Figure 4-7 shows this data superimposed on the quadrangle 
sheet. Four groups of pixels to be tested for qualification are 
circled and identified, A through D. Neither map suggest impound- 
ments on the Skagit River. Clear Lake and Big Lake are natural 
lakes, shown both on the topographic map and on the computer map. 

Target A — This coincides with Beaver Lake on the map. The 
square corners and straight-line southwest edge of this 
lake suggest that even if it was natural, it seems to 
have been modified, perhaps in the construction of the 
road shown on the map. It is also pertinent to note 
that this lake, like Clear Lake just to the north, sits 
on the same contour plane which extends from the Skagit 
all the way south to a levee. Because the contour inter- 
val on 15-minute quadrangle maps is 50 feet, there is 
the possibility of an impoundment structure of interest. 

The map does show a draining creek. A final data point 
is a cross-check in the existing National Dam Inventory. 

There is a Beaver Lake Dam in Washington, but not in 
this location. Target A needs further investigation. 

Target B — Whether these three pixels are contiguous or not 

cannot be determined from the 1:62 500 scale computer map. 

The 1:24 000 computer map in figure 4-8 shows they are con- 
tiguous . And the topographic map shows that they correspond 
to a small natural water body, Mud Lake. 

Target C — The water pixel count in this large flat area 

located south of a bluff 300 feet above the Skagit is 
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Figure 4-7.— DAM package computer-generated map 
of figure 4-6 superimposed on portion of Clear 
Lake 15-minute USGS quadrangle map shown in 
figure 4-5. 
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Figure 4-8.— Portion of DAM package computer-generated 
1:24 000-scale map corresponding to pertinent area 
of 1:62 500-scale map depicted in figure 4-7. 
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95. This represents a water surface of the order of 
150 acres, and is not indicated at all on the quadrangle 
sheet. The marginal data for the quadrangle sheet states 
that it was based on planetable surveys in 1941, updated 
by aerial photographs in 1947-48, and field checked in 
1956. The satellite data is, of course, 13 years newer. 

The elevation contours show a watershed quite capable of 
feeding an impoundment. The area is flat, and the Sun 
angle is higher than the potential shadowing of the Cultus 
Mountain system to the southeast. (If this conclusion is 
challenged, it is easy to construct a terrain profile along 
the Sun azimuth line from the top of Cultus Mountain, 
elevation 3950 feet, to the center of the suspected 
impoundment, about 430 feet, over a distance in excess 
of 3 miles. The suspect area is not shaded.) The analyst 
will report this finding as a new impoundment. 

Target D — These 18 water pixels are centered on a map feature 
which indicates an intermittent creek (light blue dashed 
line along the V-points of the contour lines) . The contours 
in the area of the suspect impoundment are not too helpful 
and again, they are at 50-foot intervals. It can be seen 
by examining any well defined stream line that the contour 
V's point upstream. There is sufficient watershed to feed 
an impoundment here, which drains to the northwest into 
Walker Creek. Without any information to the contrary, this 
target will be reported as a new impoundment. The most 
probable location of the new dam center can be estimated 
from the location of the edge pixels on the northern 
(downstream) region of the detected surface water and 
reported in the desired coordinate system (s) . 

Pixels between Big Lake and Target D — Three water pixels were 
detected to the northeast of Big Lake. The southernmost 
one appears to be on a plateau with no obvious drainage 
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pattern, in or out. If it is not a false alarm, it is 
probably a small pond. The same conclusion cannot be 
drawn about the two northern pixels since they are located 
on a natural drainage line. The symbols for these two 
pixels are located in diagonally adjacent print characters. 
Whether the pixels themselves are contiguous or not can 
only be determined by referring to the 1:24 000 computer 
map for this area. This information on water pixels, 
together with an analysis of the elevation contour lines 
and drainage from the topographic map, must then be used 
in determining the number of impoundments present and 
whether any of them warrant further checking. A final 
determination can only be made by field work or (perhaps) 
the use of aerial photography. 


4.7 REPORTING CRITERIA 

The analyst serves as a filter whenever he eliminates surface 
water which is obviously not impounded. All of the Skagit River 
contained in the Clear Lake quadrangle (fig. 4-5) is a case in 
point. Similarly, natural lakes are of no interest for the NPID 
program. With the possible exception of Beaver Lake, all of the 
lakes in the Clear Lake quadrangle are natural. None have either 
of the conditions which suggest an impounding dam; (1) a pro- 
nounced change in terrain elevation along one edge of the water 
body, determined by analysis of the map contours all around the 
water, and (2) a man-made shaped edge. The latter is readily 
discernable in large scale maps. Figure 4-9, a section of the 
Huntsville (Texas) 7.5-minute quadrangle, offers several excellent 
examples. Figure 4-10 shows the same map with computer-generated 
map data superimposed. The map was printed in 1963, based on 
aerial photographs taken in 1960, with field checks in 1963. The 
Landsat data was collected on August 29, 1972. Club Lake (1940), 
Elkins Lake (1920) , and Sunset Lake (1951) dams are all on 
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Figure 4-9.— Portion of Huntsville (Texas) 
7.5-minute USGS quadrangle map. 
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Figure 4-10.— Portion of Huntsville (Texas) 
7. 5-minute USGS quadrangle map with corres 
ponding DAM package computer-generated map 
superimposed . 
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the map and confirmed by the satellite data. Baldwin Lake, 
north of Elkins Lake, seems to have disappeared. On the other 
hand, it measures only about 150 meters wide in the Landsat 
scan direction. It is possible that no pure water pixels were 
seen by the Landsat scanner. Contours on the downstream side of 
the dam suggest that it is 20 feet or more high. This dam was 
not included in the 1974 National Dam Inventory for Walker County, 
Texas. Because the Baldwin Lake impoundment is marginal from 
both detection and qualification considerations, it should be 
reported to supervision for a subsequent decision on whether or 
not to have it field checked. The computer-generated map also 
shows a new water impoundment, well over ten acres in surface 
area, on McDonald Creek, 300 meters north of Club Lake. This 
should be reported. 

All of the preceding text calls for the following criteria for 
reporting information deduced from analysis of both conventional 
maps and computer-generated maps . 

a. Rejected pixels, ruled as not being water, should never be 
reported. 

b. All non- impounded water is excluded from reports 

c. Water impoundments which appear in the National Dam Inventory, 
or the existing map, or both, are reported. 

d* Qualifying water impoundments which appear on the published 
map, but not on the computer-generated map are reported. 

e. Qualifying water bodies which appear on the computer-generated 
map, but not on the published map are reported. 

f. Any substantial changes in the areal extent of water bodies 
should be reported. 
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4.8 REPORT ELEMENTS 

The formats and techniques for consolidating data on dams exist 
in all Corps of Engineers district and regional offices. It is 
not the purpose of this manual to codify these existing procedures. 
The report elements suggested below summarize the information 
which should develop during the correlation and interpretation 
of the computer-generated maps and published maps. 

a. Name and scale of published map 

b. Map origin coordinates (latitude- longitude of the southeast 
corner of the map) 

c. Date of last map revision 

d. Date and scene number of Landsat scene 

e. Name of dam (if known) 

f. Name of river or stream 

g. Center of interpreted impounding structure (latitude- 
longitude) 

h. Number of water pixels and/or surface area of impoundment. 

i. Is the dam shown on the published map? 

j. Is water seen on the published map missing in the computer- 
generated map? 

k. Is the dam listed in the district (regional, national) 
inventory? 

l. Area change: the impoundment shows on both the published 

maps and the computer-generated maps, but appears to be 
significantly different in areal extent and/or shape. 

4.9 INTERPRETING WITH INADEQUATE PUBLISHED MAPS 

There are regions of the United States where only 1:250 000 topo- 
graphic maps are available. Since these are unsuitable for local 
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terrain analysis, other map sources must be sought. The best 
substitute will be up-to-date orthophoto maps. If the regional 
processing center has been provided advance information on the 
scales and formats of these published maps, then corresponding 
classification maps (as well as their constituent unpublished 
1:24 000 7.5-minute maps) will be generated by the DAM software. 

4.10 INTERPRETATION WITH NO PUBLISHED MAPS 

Under these circumstances, the interpreter can only report all 
water bodies detected that do not appear to be large natural 
water bodies. Because the computer map does not convey small 
shapes well, there is no way to distinguish between impound- 
ments and small natural lakes. The analyst can only report 
locations of suspect impoundments for checking by other means. 
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